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^ , ABSTRACT 

Aims. We present multiwavelength observations of one of the most distant gamma-ray bursts detected so far, GRB 080913. Based on 

, these observations, we consider whether it could be classified as a short-duration GRB and discuss the implications for the progenitor 

^ . nature and energy extraction mechanisms. 

' Methods. Multiwavelength X-ray, near IR and millimetre observations were made between 20.7 hours and ~16.8 days after the event. 

' Results. Whereas a very faint afterglow was seen at the 3.5m CAHA telescope in the nIR, the X-ray afterglow was clearly detected 

I , in both Swift and XMM - Newton observations. An upper limit is reported in the mm range. We have modeled the data assuming a 

. collimated > 3° blast wave with an energy injection at ~ 0.5 days carrying 5 ~ 10^^ erg or approximately 12 times the initial energy 

• ' of the blast wave. We find that GRB 080913 shares many of the gamma- ray diagnostics with the more recent burst GRB 090423 for 

' being classified as short had they ocurred at low redshift. If the progenitor were a compact binary merger, it is likely composed by a 

, NS and BH. The Blandford-Znajek (BZ) mechanism is the preferred one to extract energy from the central, maximally-rotating BH. 

OO , Both the magnetic field close to the event horizon (S) and the BH mass (Mj,,,) are restricted within a relatively narrow range, such that 

' (fi/3 X 10'^ G)(Mbh/7 Mq) ~ 1. Similar constraints on the central BH hold for coUapsar-like progenitor systems if the BZ-mechanism 

^ ' works for the system at hand. 

' Key words, gamma rays: bursts - techniques: photometric - cosmology: observations 

1 . Introduction However, de Ugarte Postigo et al. (2006) observed GRB 060121 

, , „ , , • • (T9o=2s; Arimoto et al. 2006) and provided a most probable 

Gamma-ray bursts (GRBs) are generally classified mto two main hotometric redshift of z = 4.6. They suggested that this burst 

groups (Kouvehotou et al. 1993), those with short duration and ^^^j^ ^j^^ ^^.^^ ^ ^j^^^ ^^^^^ gamma-ray bursts residing at 

hard spectra and those with ong duration and soft spectra, j^^gj^ ^^^^^^^^^ ^j^^^j^ p^^^^^j^ ^^j^^^^^ 

to a different progenitor 

This simplistic classification scheme could be more complex, as ^^^^^ high-redshift bursts were later studied statistically 

shown by several studies (Zhang et al. 2009; Horvath et al. 2006, g^^g^^ ^1 (2007b), who found that a significant number of 

. . . distant short bursts could exist. 

For short bursts the general idea is that they originate in the 

near Universe, at redshifts significantly lower (z~0. 5) than those „„„ nonr,^^ j- j u c n c onno 

f 1 /-Dt> / r-DD nfnsmu . a oot di . i 080913 was discovered by Swift on 13 Sep 2008 

of long GRBs (e.g. GRB 060502b at z = 0.287, Bloom et al. ,c u a » i oaaon -tu u » ^ ^ a ^ ttt- j 

onm %r,v> ncmi . n ^a^a c a u /i oaa^; /-dd (Schady et al. 2008). The burst started at 06:46:54 UT and 

2007; GRB 051221a at z = 0.5464, Soderberg et al. 2006; GRB , , o i • . ^ . • • .u < j 

, „ 1^.^ n , ^^^^ >,r.Ti r^cr^-^r^A lasted for 8 s, placing it, at nrst sight, in the lone-dumtion 

050911 at z = 0.1646, Berger et al. 2007a; GRB 050724 at z , f nou /c. .i \ i onncf n i u au 

rv-icT Ti . 1 .-irvrvc /-^nn ncnnnf^ . rv i ^rv t- class of GRBs (StamatiKos ct al. 2008). It was also obscivcd by 

= 0.257, Berger et al. 2005; GRB 050709 at z = 0.160, Fox „ mn\Tr^ au a a f , n i a-7 -2 

» 1 onnc Av>T, n'^f^'^rKau , a n-it. /- u i . i oAns:\ KonusfWlND and had a fluence of (5.6 + 0.6) x 10 erg cm ^ 

et al. 2005; GRB 050509b at z = 0.226, Gehrels et al. 2005). ■ ^. ic icm i • r^nn t .u 

in the 15-150 keV range, making It an average GRB. In the com- 

offprint requests to: D. Perez-Rami'rez, e-mail: bined BAT- WIND spectrum, the observed prompt energy spec- 

dperez@ujaen.es trum could be best-fitted by a power-law with an exponential 

* Based on observations taken with the 3.5m Calar Alto telescope, cutoff model dN/dE ~ X exp(-(2H-a) X E/Ep,ak) with a = 

the Plateau de Bure interferometer and the XMM - Newton satellite. -0.89;^Q 4g and Ep^ak - 131^48^ keV (Pal'shin et al. 2008). 
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The prompt dissemination (21 s) of the GRB position by 
Swift enabled instant responses of robotic telescopes, such as 
the REM robotic telescope (D'Avanzo et al. 2008). Rapid ob- 
servations obtained by GROND at the 2.2m telescope in La Silla 
allowed the identification of a near-infrared (nIR) counterpart 
(Rossi et al. 2008) ~ 3 min after the burst trigger. Shortly after 
Swift slewed and started data acquisition, a fading X-ray source 
was detected by the ^wZ/f/XRT, which was identified as the GRB 
080913 afterglow (Beardmore et al. 2008). This triggered a mul- 
tiwavelength campaign at different observatories aimed at study- 
ing the afterglow. A photometric redshift in the range 6.1-6.7 
was derived (Greiner, Kruehler and Rossi 2008), and a spectro- 
scopic z - 6.7 was later confirmed by a VLT spectrum (Fynbo 
et al. 2008, Greiner et al. 2009). This implied rest frame val- 
ues of T90 ~ 1 s and Epeakjesr ~1000 keV consistent with a 
short- duration GRB (Pal'shin et al. 2008). This is also supported 
by the negligible spectral time lags found in the BAT energy 
range (Xu 2008). For a standard cosmology model with Ho - 
71 km/s/Mpc, Qm = 0.27, Qa = 0.73, the isotropic energy re- 
lease is Eiso ~ 1 X 10^^ erg (1 keV-10 MeV, rest frame), with a 
look-back time of 13.67 Gyr 

A more recently detected burst, GRB 090423 (Tanvir et al. 
2009, Salvaten-a et al. 2009), is a further extreme redshift (z=8.2) 
and potentially short GRB which exhibits similar properties, 
such as burst duration (Tgo ~ 10.3 ± 1.1s, and a rest-frame dura- 
tion of ~ Is), spectral lag times negligibly small, consistent with 
zero, and an intrisincally hard spectrum as GRB 080913. 

This burst, together with GRB 080913, points to the fact that 
the current dichotomy is not always consistent. That is, standard 
indicators of the physical nature of GRBs, such as duration and 
hardness, may no longer be the only diagnostic used in phys- 
ically classifying high-redshift GRBs. They reveal a need for 
a revision of the traditional observational criteria. Zhang et al. 
(2009) tackled this question proposing new operational proce- 
dures in the determination of the physical category of GRBs. 
According to this work, GRB 080913 and GRB 090423 are con- 
sidered to belong to the Type II category (i.e. massive-star core 
collapse origin). However, these two high-z bursts may also 
be compatible with a "specific Type I scenario" driven by the 
Blanford-Znajek mechanism in a BH-NS merger. 

Here we report multiwavelength observations carried out, 
from the millimetre to the X-ray band, in order to study the 
afterglow of GRB 080913. We also discuss the implications of 
these observations for the nature of short-duration GRBs progen- 
itors. Finally, we include some results extracted from the litera- 
ture about the implications for the nature of long-duration GRBs 
progenitors 

2. Observations and data reduction 

2.1. Near-IR observations 

Target of opportunity (ToO) observations in the nIR were trig- 
gered starting 20.7 hr after the event with the 3.5 m telescope 
(-1-OMEGA 2000; Bailer-Jones et al. 2000) at the German- 
Spanish Calar Alto Observatory (CAHA). A 5,500 sec image 
was acquired in the /-band filter (see Fig. 1), with a 1.35" aver- 
age seeing. We followed the standard data reduction procedures 
such as dark and sky frame subtraction, plus bad-pixel mask and 
master flat-field correction. The photometry for our final image 
was performed by means of the PHOT routine under IRAfQ A 



' IRAF is distributed by the NOAO, which are operated by USRA, 
under cooperative agreement with the US NSF. 




Fig. 1. The 7-band image of the GRB 080913 field taken at the 
3.5m CAHA telescope on 14 Sep 2008 (20.7 hr after the burst 
onset). The error circle marks the XRT position of the counter- 
part with a radius of 1.9" (Beardmore et al. 2008). 



range of apertures were checked, and the one yielding the mini- 
mum photometric error was selected. The candidate initially re- 
ported by Rossi et al. (2008) was barely detected (2.3o- level) 
in the /-band image (Fig. 1), with an estimated J-band Vega 
magnitude of 22.4+0.5, including the calibration zero point error 
(0.15 mag) given by the 2MASS Catalogue. In order to compare 
our detection with the lightcurve for this GRB as presented by 
Greiner et al. (2009), we have evaluated the AB-to-Vega system 
magnitude offset coefficient for the OMEGA 2000 instrument 
in the /-band to be 0.97. The converted J-band magnitude in the 
AB system for our detection is 23.4+0.5, which at the time of the 
observation, agrees with the GROND data (Greiner et al. 2009). 
Our detection occurred close to the peak of the re-brightening 
phase (Fig. nil. 

2.2. Millimetre observations 

Additional mm observations were obtained at the Plateau de 
Bure Interferometer (PdBI) as part of our ToO programme. 
The PdBI observed the source on different occasions in the pe- 
riod of time of three days in compact configuration. We used 
the carbon star MWC349 as primary flux calibrator (assuming 
F(v) = 1.1 ■ (v/86.2GHz)" '') with the amplitude and phase caH- 
brations relative to the quasar 0454-234. The data reduction was 
done with the CLIC and MAPPING software distributed by the 
Grenoble GILDAS group. We analysed the data with position- 
fixed (RA(2000) = 04:22:54.66, Dec(2000) = -25:07:46.2) fits 
in the UV plane, which only yielded upper limits but allows to 
constrain the mm-lightcurve. The 3cr-limits are 0.72 mJy (99 
GHz, Sep 16.1 UT), 1.44 Jy (84 GHz, Sep 21.2 UT) and 0.90 
mJy (106 GHz, Sep 30.1 UT). 

2.3. X-ray observations 

We made use of public Swift/XRT data obtained starting ~94s 
after the BAT trigger time (Tq). The data were acquired in the 
XRT Photon Counting (PC) mode. 

The lightcurve in the 0.3-10 keV energy range (Evans et al. 
2007) exhibits a number of flares in the initial orbit, with the 
peak of the largest flare observed at To+1800 s (in the observer 
frame). The decay in the time interval To +400 s to To+llOO s 
can be fitted by a power-law with an estimated decay index ax 
- l-2;^Qj (where fx oc r""). The spectrum corresponding to the 
period of flare activity in the initial orbit, can be fitted by an 
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Fig. 2. The EPIC-pn image of the high redshift GRB 0809 1 3 ob- 
tained with XMM-Newton ~4.5 days after the burst onset. The 
image has been adaptively smoothed to emphasize the detec- 
tion of the GRB. The position derived by Swift for this burst 
is marked by the circle. 



absorbed power-law with a photon index of F ~ 1-7^" 4, assum- 
ing that the Galactic column density value in the direction of the 
burst is 3.2 x 10^" cm^^ (Kalberla et al. 2005). 

A ToO XMM-Newton observation (Obs. ID. 0560191701) 
started on Sep 17.61 UT, i.e. 4.3 days after the 5w//r/BAT trig- 
ger The European Photon Imaging Camera (EPIC) CCD cam- 
eras were operated in the Full Frame Mode. The total exposure 
time for the EPIC-pn camera was 14.0 ks. The Thinl optical 
blocking filter was used for the EPIC-pn camera, whereas the 
medium filter was used for both EPIC-MOS cameras. The XMM- 
Newton Observation Data Files (ODE) were processed using 
XMM-Newton Science Analysis Software (SAS version 7.1.0) 
and the calibration files from the Calibration Access Layer as on 
14 Dec 2007. After excising periods of high-background, the net 
exposure times of the EPIC-pn, MOSl, and MOS2 observations 
are reduced to 3.6 ks, 10.8 ks, and 7.9 ks, respectively. 

The XMM-Newton EPIC observations of GRB 080913 de- 
tected X-ray emission from its afterglow (see Fig. 2) at EPIC- 
pn, MOSl, and MOS2 count rates of 0.0109+0.0021 cnts s-\ 
0.0032+0.0008 cnts s"', and 0.0019+0.0008 cnts s"', respec- 
tively. The net count number, ~90 counts, is not sufficient 
to carry out a spectral fit. Adopting an absorbed power-law 
model of spectral index F=1.7 and column density A^h = 3.2 x 
10^*^ cm"^, the EPIC spectral shapes and count rates imply an ab- 
sorbed X-ray flux of 3.8x10"''* ergs cm"^ s"' and an unabsorbed 
X-ray flux of 4.1x10"''* ergs cm"^ s"' in the energy band 0.3-10 
keV. The X-ray luminosity in this same band at restframe is 5.7x 
lO''^ erg s"' (assuming in XSPEC Hq = 71 km/s/Mpc and Qa = 
0.73). 



3. Results and discussion 

Following the discovery of the X-ray afterglow with Swift, we 
detected a faint nIR afterglow (consistent with the position given 
by Rossi et al. 2008) and the X-ray afterglow 4.5 days after with 
XMM-Newton. We discuss in this section the classification and 
likely progenitor of this burst. 
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3.1. Spectral Flux Distribution of GRB 080913 

Using the model and methods described by Johannesson et al. 
(2006) we fitted the multiband observations of the afterglow to a 
fireball model with energy injections. 

One injection is needed in order to account for the bump seen 
at ~ 0.5 days in the light curves. Since no jet-break is seen in the 
light curves up to around 10 days, we can only put a lower limit 
to the collimation angle of the jet, 60 > 3°. This is delayed from 
the definition of Sari et al. (1999) due to the energy injection, 
the sideways expansion of the jet and the detailed calculation of 
the equal arrival time surface. Our preferred scenario (giving the 
best fit) for a collimation of 3° is an initial energy release of £0 
= 4x 10^' erg into a uniform medium with density no = 2 cm"-'. 
This fit results in a;^'^/d.o.f = 170/50 where the high value is 
mainly caused by the scatter in the X-ray light curve the model 
is unable to reproduce. In order to explain the bump at 0.5 days, 
an energy injection carrying approximately 12£'o (5x 10^^ erg) is 
needed. Due to the lack of a turnover in the light curves at early 
time, we require the initial Lorentz factor of the blast wave to be 
Fo > 500. The SFD is best fit with an electron index p -2.11 and 
micro- physical parameters e, - 6x10"'' and eg - 3x10"^. Note 
that the parametrization of the electron population has changed 
from Johannesson et al. (2006) and we now follow the tradition 
of Panaitescu & Kumar (2001). The minimum Lorentz factor of 
the electron distribution is now defined as ymin - e;7«p/?Me(r- 1). 
Figure [3] shows the radio to X-ray SFD predicted by our model 
for 3 epochs together with observation data corrected for intrin- 
sic extinction. The Galaxy extinction is negligible. 

Please note the the above values for the best fit parameters 
are in many cases highly dependent on the value chosen for 6q. 
The energy required for the blast wave goes as 0^ and lower val- 
ues of 60 put an upper limit on «o from the jet break require- 
ments. Additionally, there is a strong correlation between the 
value of «() and eg because we do not have the required data 
to constrain the SFD at lower frequencies. Other parameters are 
less sensitive to the value of ^o, especially p that is well con- 
strained from the nIR SFD. 



3.2. Comparison with high-redshift long GRB 050904 and 
short GRB 060121 

Given their highest redshifts, one could quaUtatively compare 
the properties of GRB 080913 (z = 6.7), GRB 060121 (z = 4.6) 
and GRB 050904 (z = 6.295). The reason for choosing the latest 
cases is that GRB 060121 is a short duration GRB (T9o~2s), very 
likely at high-redshift (z - 4.6, de Ugarte Postigo et al. 2006), 
whereas GRB 050904 is a long GRB (T9o~3 1 s ) at a comparable 
redshift (z = 6.295; Haislip et al. 2006; Kawai et al. 2006). We 
have constructed the restframe isotropic 0.3-10 keV luminosity 
lightcurves of the three GRBs, assuming a power law spectrum 
with a photon index F. The time evolution of F was found by 
using linear interpolation between the F values determined from 
Swift /XRT spectra on a logarithmic time scale. This was done to 
get a smooth K-correction with time. As seen in Fig.|4]the 0.3-10 
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Fig. 3. Spectral Flux Distribution of the afterglow from radio to 
X-rays 0.01, 0.10, and 3.09 days after the burst in the observer 
frame. Filled circles are based on the nIR data from Greiner et al. 
(2009) and the X-ray data from SwiftpiKY. Triangles represent 
upper limit: in the Tq-hO.OI days plot they correspond to nIR 
observation upper limits whereas in the To+3.09 days plot it is 
the 3-cr upper limit obtained at PdB (this paper). 
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3.3. Implications for a short duration progenitor origin 

Because of the intrinsically short duration, hard spectrum and the 
effectively zero spectral lags of GRB 080913, we consider in the 
following the constraints on the typical progenitor systems in- 
voked for standard short-hard GRBs. We note that a late energy 
injection, used to fit the afterglow light curve (Sect. 13.1b . does 
not necessarily argues against a short GRB progenitor or, equiv- 
alently, imply a long lived central engine, if the ejected material 
has a distribution in initial Lorentz factors. Under the assumption 
that the initial material powering the afterglow was ejected with 
a Lorentz factor of ~ 500, matter ejected simultaneously with 
a Lorentz factor of ~ 40 would refresh the afterglow at around 
0.5 days. The rather large energy injection of 10 times the initial 
energy release requires, however, a significant increase in the 
efficiency of the central engine that poses constraints on either 
short or long lived central engines. 

At a redshift z = 6.7, the burst ocurred when the Universe was 
only ~ 0.8 Gyr old. According to Yoshida et al. (2008), the first 
stars formed ~ 0.2 Gyr after the Big Bang. Thus, if GRB 0809 1 3 
resulted from a merger of compact objects, the whole evolution 
of the progenitor stars (Tevoi) and the merger time (Tmer) ha- 
ppened within T = Tevoi + T"mer ^ 0.6 Gyr Such a small value 
of T sets a lower limit to the mass of the progenitor stars. Thus, 
neither of the two stars forming the binary system could be less 
massive than ~ 4Mq. This only restricts the type of progenitor 
system if the burst were a result of the merger of a white dwarf 
(WD) and a BH0 

As the lower accretion timescales (~ 10 s) occur for the 
higher-mass WD (Popham, Woosley & Fryer 1999), the WD 
should have a mass close to the Chandrasekar mass (~ IAMq). 
For NSh-NS/BH binaries, the evolution of the progenitor stars is 
only a small fraction of r (Tgvoi ~ 10-20 Myr; Belczynski et al. 
2006). Thus, T roughly equals Tmer in this case. 

Since Tmer ~ 0.6 Gyr, if the GRB originated from a double 
NS (DNS) merger, it could have proceeded through a classical 
merger channel, for which Tmer ~ 0.1 - 15 Gyr (Belczynski et 
al. 2006). Considering the large of GRB 080913, it is more 
likely that it was hosted in a NS-hBH merger than in a DNS. The 
latter type of mergers have typically t^I^^^ ~ 1 Gyr, but with a 
non-negligible merger probability for times O.lGyr < tJ^^^^^^ < 
IGyr 

According to Oechslin & Janka (2006) the total equivalent 
isotropic energy released in gamma-rays (£,so) from a progen- 
itor consisting of a BH-accretion torus system, which boosts a 
neutrino driven ultrarleativistic jet is 



(1) 



Fig. 4. The K-corrected (rest frame) 0.3-10 keV isotropic lumi- 
nosity evolution of GRB 080913 (red) in comparison to the long- 
duration GRB 050904 (blue, z - 6.3) and the short-duration 
GRB 060121 (green, z - 4.6). The solid square at f - Icrb - 
4.9 X lQ*s represents the XMM-Newton observations carried out 
for GRB 080913. As seen GRB 080913 shows lower isotropic 
luminosities than the short GRB 060121. 



keV decay of GRB 0809 1 3 is similar in character to those exhib- 
ited by GRB 060121 and GRB 050904. However, its isotropic 
luminosity at early stages is lower by a factor ~30. Therefore, 
based on the afterglow lightcurve it is not possible to rule out 
any of the possible origins for the progenitor. 



where Mace and c are the accreted mass and the light speed, 
respectively, /i, /a, /a, and /4 are different efficiency factors 
(see below), and /n - O^^ll is the jet collimation factor For 
GRB 080913, both fa and f^, the fraction of the energy of ul- 
trarelativistic jet matter, which can be emitted in gamma-rays in 
course of dissipative processes that occur in shocks, can be fixed 
from the light curve fit shown in Sect. 13. II A jet collimation fac- 
tor /n = 1.37 X 10"^ results taking as jet half-opening angle 
= 3° (Sect. 13.1b . This value is smaller than, but consistent 
with /n ~ 0.015 - 0.034 obtained in numerical models (Aloy, 
Janka & Miiller 2005), since 6q is only a lower bound in our 
case. The factor /4 is set by the ratio - fnEy/EAc = 0.024, 
where Eac is taken here equal to the kinetic energy used to 



- WD-l-NS mergers possess much longer timescales and are more 
likely progenitors of long GRBs (e.g., King, Olsson, & Davies 2007) 
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model initial afterglow in Sect. 13.11 i.e., Eag - E^. Such fig- 
ure is also consistent with estimates of the internal shock model, 
/4 < 0.3 (e.g., Mimica & Aloy 2009 and references therein). For 
the remaining factors in Eq. ([TJ, we find that the observed en- 
ergy in the prompt GRB phase needs the concurrence of a large 
accretion disk mass {Mace ^ O.SSMq; Oechslin, Janka & Marek 
(2007) obtain Mgcc ^ 0.3Mo) and several large efficiency fac- 
tors (probably so large that they rule out a neutrino mediated 
energy extraction from the central engine). A large conversion 
efficiency of the accreted mass into neutrino emission f\ > 0.1 
(to be compared with a typical value f\ > 0.05; e.g. Setiawan 
et al. 2006; Lee et al. 2005), a conversion efficiency of neutri- 
nos and antineutrinos by annihilation to e"^ pairs /2 > 0.06 (for 
reference, /2 ~ 0.001, . . . , 0.04 is estimated by Ruffert & Janka 
(1999), Setiawan et al. (2006), or Birkl et al. (2007), and a large 
fraction of the -photon fireball energy which drives the ultra- 
relativistic outflow > 0.4 (larger than ~ 0.1 extracted from 
simulations; Aloy, Janka & Miiller 2005). The combination of a 
large /2 and a large disk mass is supported by the steady models 
of Birkl et al. (2007) with values of the dimensionless angular 
momentum of the central BH a ~ 0.4 - 0.5. 

We note that if Eag is taken to be the sum of the contribu- 
tions due to the initial kinetic energy (£0) and of the late energy 
injection (12£'o; Sect. 13.1b . then /4 ~ 1.3 x 10"^, implying that 
all the remaining /ree efficiency factors (/i, /a and /a) have to 
be, at least, a factor of 5 larger than the largest estimates of them, 
which suggests that a neutrino driven outflow cannot account for 
the observed phenomenolgy. 

We point out that the previous analysis is only sensitive to the 
value of 6() indirectly, if we seek to accommodate with our model 
the values of Eag- This is because Eag - fififiMaccC^, but 
the fitted values of Eag depend on 0^ (Sect. [TTJ. Berger (2007) 
makes the hypothesis that outflows of short-duration GRBs with 
the highest energies are strongly collimated. If the outflow is 
neutrino-driven, we point out that such hypothesis seems consis- 
tent with the fact that the largest Eis,, are linked to BH-torus sys- 
tems in which either the torus mass is large or /i is large. A large 
torus mass may arise in mergers between compact objects of 
different masses (Oechslin & Janka 2006; Shibata & Taniguchi 
2008). The conversion efficiency of the accreted mass into neu- 
trino emission tends to increase with increasing values of the 
torus viscosity (Setiawan, Ruffert & Janka, 2006). Remarkably, 
a large viscosity yields more vertically extended (inflated) ac- 
cretion tori, which may help to collimate the ultra-relativistic 
outflow in narrow channels. 

Given the large efficiencies needed to account for the large 
Eiso of GRB 080913 if the outflow were neutrino driven, it seems 
more natural in this subclass of extremely energetic short GRBs 
to invoke an energy extraction mechanism directly linked to the 
BH spin, e.g., the Blandford-Znajek (BZ) process (Blandford & 
Znajek 1977). In this case, an estimate of the total power pro- 
duced by the central engine is (Lee, Wijers & Brown 2000), 
Pbz = 1.7 X lO^'^a\B /lO^^Gf (Mbh/Maffih), wher e f(h) ^ 
[(1 + h^)/h^][(h + l//z)arctan(/i) - 1], // = a/(l + yjil -a^), 
and B is the magnetic field strength at the event horizon of the 
BH with mass M/,/,. We obtain the isotropic equivalent energy 
released in y-rays due to this process during the intrinsic event 
duration (Tgo ^ 1 s) as E/so = Pbz * 7^90 * /s * /4- Note that in the 
previous estimate, we use the same efficiency factors ~ 0.4 
and /4 = 0.024 as in the previous paragraphs to account for the 
facts that (1) only a fraction of the released energy will be used 
to drive a ultrarelativistic outflow, and that (2) the radiated en- 
ergy in y-rays is much smaller than the kinetic energy of the 



outflow Eag- Hence, the observed energy can be reached if ei- 
ther the value of the dimensionless angular momentum of the 
central BH is a ~ 1, the magnetic field surrounding the BH is 
B > lO^'G or the BH has a mass Mbi, > 20Mq. We note that 
to form a 20 M© BH in a low-metallicity star (which shall be the 
case at the redshift of GRB 080913), the initial mass of such star 
shall be > 50 Mq (Woosley, Heger & Weaver 2002). 

Since Pbz depends quadratically on both B and M/,/,, we 
may estimate which is the range of variation of these two pa- 
rameters such that the resulting Ejso complies with the energet- 
ics observed for GRB 080913. Lower values of the BH mass, 
Mhh ~ 3Mq, closer to the typically considered ones in mergers of 
compact objects, require extremely large values of the magnetic 
field strength (B > 7 x lO'* G). Such magnetic fields would prob- 
ably brake excesively the rotation of the stellar progenitor core, 
likely inhibiting the formation of a maximally rotating Kerr BH 
(i.e., reducing the value of a sensitively below 1). Hence, in ac- 
count of the large estimated mass of the BH, if a merger of com- 
pact objects were the progenitor system of this GRB, a NSh-BH 
merger is favoured, since the typical mass of the BH resulting 
from a DNS merger is < 3Mq. On the other hand, if we consider 
values of the magnetic field strengh smaller than the reference 
value of 10'^ G, we find that even a factor of three smaller field 
yields a BH mass Mj/, ~ 60 M©. To form such a massive BH the 
progenitor star should have a mass ~ I4OM0, i.e., in the limit 
of being stars which are disrupted at the end of their lives by 
the pair instability without leaving any remnant BH (Woosley, 
Heger & Weaver 2002). Thus, both the magnetic field strength 
and the BH mass are restricted within a relatively narrow range, 
3 X 10'^ G < B < 3 X lO'*" G, and 7 Mo < Mm < 60 Mq, respec- 
tively. Within the former range of values for the triad of param- 
eters a, B and M/,i, the most likely ones are those favouring the 
largest possible value of a. Because of the fact that strong mag- 
netic fields tend to slow down the rotation of the stellar core, and 
because of the difficulty to build up magnetic fields in excess of 
~ 10'^ G by the collapse of stellar cores (e.g., Obergaulinger, 
Aloy & Miiller 2006, Obergaulinger et al. 2006), even consider- 
ing the action of the magnetorotational instability (Akiyama et 
al. 2003; Obergaulinger et al. 2009), values of a ~ 1 fit bet- 
ter with the lowest values of B in the aforementioned range. 
Thereby, to reach the appropriate BZ power to fuel GRB 0809 13, 
we favour BH masses in the upper end of the interval stated 
above. 

As noted above, if the total energy of the afterglow is Eag - 
l3Eo, the reduced value of /4, yields even more stringent con- 
straints on the central engine, since Pbz has to be 13 times larger, 
which needs of a ~ I, B > 1.6 x 10"'G and Mw, > 50Mo. 
However, we have to be cautious with the inferences based on 
the values of the total kinetic energy in the afterglow, since dif- 
ferently from Ejso (directly measured), Eag results from a model 
fit of the afterglow light curve, which sensitively depends on the 
value of 00 (SectlTTTi. 

FinaUy, the large redshift of GRB080913 fits in theoretical 
models where the rate of NS+NS/BH mergers follows either the 
star formation rate or the star formation rate with delays smaller 
than 1 Gyr (Janka et al. 2006). 

3.4. Implications for a long duration progenitor origin 

Recent studies (Zhang et al. 2009, Belczynski et al. 2009) have 
suggested that GRB 080913 has a long-duration progenitor ori- 
gin. 

Zhang et al. (2009) present a new scheme for classifying 
bursts based on criteria more closely related to the progenitor 
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type to differentiate physical origins, such as SN association, 
host galaxy properties and the offset of the GRB location in the 
host galaxy. They classify bursts into two main categories: Type I 
(with a compact star-merger origin) and Type II (with a massive- 
star core collapse origin). Under such scheme, GRB 080913 and 
GRB 090423 are identified as Type II candidates based on (i) 
the geometrically-corrected gamma-ray (Ey) and kinetic (Ek) 
energies (with large values), (ii) intrinsic afterglow luminosi- 
ties (moderately bright), (iii) the high density of the circumburst 
medium and, (iv) the marginal compliance of the Ep - Eyjs,, re- 
lation. However, because these bursts are intrinsically short, but 
still considered more likely long, they suggested a possible "spe- 
cific Type I scenario" driven by the Blandford-Znajek mecha- 
nism of a BH-NS merger 

The analysis made in the previous section, regarding the 
properties of the central BH in order to be able to deliver the 
sought Eis„ of GRB 080913, is formally independent of the fact 
that the progenitor system is a single massive star or a mem- 
ber of a binary. Thus, we cannot exclude the possibility that the 
progenitor system is a massive low-metallicity star that forms a 
collapsar-like engine (or a Type II GRB according to the classifi- 
cation of Zhang et al. 2009) whose energy is extracted by means 
of a BZ-mechanim. In such a case, the estimated BH mass and 
magnetic field are the same as in Sect. 13.31 

Belczynski et al. (2009) point out that based on the currently 
used gamma-ray diagnostics (T90, Ep^ak, hardness ratio) these 
bursts would be considered as short had they occured at low red- 
shift. They argue that based on the average Swift detection rates 
(accounting for selection effects), estimated for long GRBs to 
exceed 10 times the rates for short GRBs, these bursts might 
belong to the long class. At the redshifts of these GRBs, the cal- 
culated rates are 1 yr ' and 0.1 yr ' per unit redshift, for long 
and short GRBs, respectively. 

4. Conclusions 

We report multiwavelength observations of the high z, poten- 
tially short-duration gamma-ray burst GRB 080913 acquired 
between 20.7 hours and 16.8 days after its detection by Swift. 
Although nIR and X-ray afterglows were found, no mm after- 
glow was detected. The X-ray spectrum is consistent with neg- 
ligible intrinsic absorption. We have modeled the data with a 
collimated (do > 3 deg) blast wave with an energy injection at 
0.5 days, requiring a total energy release of more than 5 x 10^^ 
erg. 

At a redshift z = 6.7, the burst ocurred when the Universe 
was only -0.8 Gyr old. If GRB 080913 resulted from a merger 
of compact objects, the whole evolution of the progenitor stars 
and the merger time happened within < 0.6 Gyr, following the 
star formation rate with delays of < 1 Gyr (Janka et al. 2006). 
This implies that neither of the two stars could be less massive 
than ~ 4Mq. If the system of compact binaries were formed by a 
WD and a BH, the short duration of the burst restricts the mass of 
the WD to be close to the Chandrasekar mass. If the progenitor 
system of this GRB is a merger of compact objects, a NS-hBH 
system is favored over a DNS mergers, having the BH a rather 
large mass (M^/, > 20 Mq) and being threaded by huge mag- 
netic fields (B > 10'^ G). An electromagnetic energy extraction 
mechanism directly linked to the BH spin, like the Blandford- 
Znajek process, may likely operate in this subclass of GRBs. 
We note that only a factor of three smaller value of B, would 
request Mhi, > 60 Mq to explain the observed Eiso- However, 
to form so massive BHs the progenitor star shall have a mass 
~ 140 Mq, i.e., they are in the limit of being stars whose final 



fate is to be disrupted by a single pulse due to the pair instabil- 
ity, without leaving any remnant BH (Woosley, Heger & Weaver 
2002), i.e., without the possibility of producing a GRB engine. 
Considering the fact that magnetic fields in excess of ~ 10'^ G 
are difficult to build up by the collapse of the core of the pro- 
genitor star, and that such fields may reduce the angular momen- 
tum of the resulting BH, our preferred choice of parameters is 
a combination of magnetic field strenght B > few x 10'^ G and 
Mhi, > 50 Mq. In regard of these figures, we temptatively suggest 
that systems with large BH masses (close to the limit set by the 
pair instability) and magnetar-magnetic field strengths may con- 
stitute the specific Type I scenario invoking the BZ-mechanism 
in a NS-i-BH mergers ocurring at high-z considered by Zhang et 
al. (2009). However, we cannot rule out the possibility that a sin- 
gle massive low-metallicity star (Type II scenario), rather than a 
NS-hBH merger, that yields a BH with the aforementioned prop- 
erties, and which fuels an ultrarelativistic ejecta by means of a 
BZ-mechanism, constitutes the central engine of GRB 080913. 

It has been also demonstrated that the synergy between mis- 
sions like Swift, automated instruments and robotic observato- 
ries, facilitate further study of the population of high-z GRBs, 
and help to assess whether a significant fraction of short-duration 
GRBs coexist at these high redshifts. 
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